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ABSTRACT 

The spectral properties from radio to optical bands are compared between 
the 18 optically bright Gamma-ray burst afterglows and well established power- 
spectrum sequence in Blazars. The comparison shows that the afterglows are 
well agreement with the well known Blazar sequence (i.e., the z/L^(5GHz)-aRo 
correlation, where «ro is the broad-band spectral slope from radio to optical 
bands). The afterglows are, however, clustered at the low luminosity end of the 
sequence, which is typically occupied by high frequency-peaked BL Lac objects. 
The correlation suggests that Gamma-ray burst afterglows share the similar emis- 
sion process with high frequency-peaked BL Lac objects. We further identify a 
deviation at a significance level larger than 2a from the sequence for three typical 
optically "dark" bursts. The deviation favors a heavy extinction in optical bands 
for the "dark" bursts. The extinction Ay is estimated to be larger than 0.5-0.6 
magnitude from the i/L^(5GHz)-aRo sequence. 

Subject headings: gamma-ray burst: general — BL Lacertae objects: general - 
methods: statistical 



1. INTRODUCTION 

With our current knowledge, Gamma-ray bursts (GRBs) and Blazars (a collective term 
for BL Lac objects and flat-spectrum radio quasars) share the similar observational phe- 
nomenon involving a relativity jet on the line-of-sight of observers. In addition to the ob- 
servational phenomena, both classes of the objects share similar central engine: accretion 
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onto a black hole with rapid spin, and similar emission processes: synchrotron and possible 
inverse Compton process (see Zhang 2007 for reviews; Urry & Padovani 1995). In spite of 
these similarities, the two classes of the objects differ from each other in following several 
respects: central engine mass, bulk Lorentz factor and power-spectral energy distribution 
sequence. 

It is now well known that the relativity jets are launched from a supper massive black 
hole with mass ~ 1O 8 ~ 9 M for Blazars according to their spectroscopic observations (e.g., 
Lacy et al. 2001). Nevertheless, the popular "collapser" model indicates that Long-GRBs 
are produced by the death of young massive stars with masses > 25M (see Woosley & 
Broom 2006 for a review). 

Direct measurements of bulk Lorentz factors (r) of jets are available for some Blazars 
through the VLBI observations of the superluminoal velocity. The observations reveal a 
typical Lorentz factor T ~ 10 for Blazars (e.g., Jorstad et al. 2001; Savolainen et al. 2010). 
However, the GRB fireball model requires a larger initial T (> 100) to ensure the 7-ray 
photons are transparent to the 77 — > e ± process (e.g., Piran 1999). 

The power-spectral sequence has been well established for Blazars by Fossati et al. 
(1998), and was recently confirmed by the Fermi-detected Blazars (Sambruna et al. 2010), 
in which the two peaks of spectral energy distribution (SED) from radio to 7-rays shift toward 
lower frequency for more luminous objects. This sequence can be quantitatively described 
by the tight correlations between the photon frequency of the first peak of the SED and 
the luminosities in different bands. Several correlations have been found for GRB prompt 
emission in recent years (see Schaefer 2007 for a summary for these correlations). Generally 
speaking, these correlations suggest that GRBs with more powerful prompt emission (or 
total released energy) tend to have higher peaked photon energy in their time averaged 
7-ray spectra (e.g., Amati et al. 2009). 

In this letter, we compare GRB's radio-optical afterglow emission with the Blazar se- 
quence. The study is motivated by the fact that the Lorentz factors T in GRBs decrease 
with both time and radius of the jets as a powerlaw, which means the values of T in the late 
time afterglows could be comparable to Blazars. The A cold dark matter cosmology with 
parameter h = 0.73, Vl m = 0.24 and = 0.76 is assumed throughout our calculations. 

2. SAMPLES 

Frail et al. (2003a) presented a complete catalog of radio afterglows of 75 pre-swift 
GRBs (between 1997 and 2001). The catalog contains the detections in the frequencies from 
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1.4 GHz to 650 GHz for 25 afterglows. In order to provide a full examination, we further 
performed a complete search from the GCN circulars and literature for the radio afterglow 
detections for the bursts triggered later than 2001. The corresponding optical light curves 
were also compiled from the public papers. Our GRB sample finally consists of the radio- 
optical afterglow detections of 18 GRBs with measured redshifts. The sampld3 is tabulated 
in Table 1. For each GRB, Column (1) and (2) lists the GRB identification and measured 
redshift, respectively. The time in days at the observer frame since the trigger is listed 
in Column (3). The next two columns present the specific fluxes at 5GHz and at 5500A 
(V-band) that were observed at the time given in Column (3). The optical V-band flux is 
corrected for the corresponding Galactic absorption, but the local one, because the statistical 
studies indicate that the GRB optical afterglows have very small averaged extinction around 
(A v ) ~ 0.2 mag (e.g., Kann et al. 2007). 

The last column presents the calculated fc-corrected monochromatic luminosity at 5 GHz. 
The monochromatic luminosity is derived from the observed flux f v listed in Column (4) as 
z/L l/ (5GHz) = \-Kd\vf v {\ + z)~ a ~ x , where di is the luminosity distance and a the spectral 
slope defined as f v oc v~ a . A typical slope a = 1/3 is adopted in the fc-correction for the 
slow cooling, optically thin phase (Sari et al. 1998). The calculated uL u (5 GHz) distribute 
in a very narrow range roughly from 10 40 to 10 42 ergs s _1 , which quality agrees with the 
similar distribution of the luminosities at 8.5 GHz (Frail 2003). 

Fossati et al. (1998) compiled a sample of Blazars with well sampled SEDs from radio 
to 7-ray. There are totally 110 Blazars with measured redshifts. These objects are used as 
a comparison sample in our subsequent analysis. The authors fc-corrected all the reported 
specific fluxes to the rest frame according to the empirical spectral slopes. Using the fc- 
corrected fluxes, we then calculate the corresponding luminosities at 5 GHz for the Blazars 
by the same method mentioned above. 



3. RESULTS 

Using the calculated average radio-to-7-ray SEDs, Fossati et al. (1998) identified a 
Blazar sequence in which there is a systematic trend of the SEDs as a function of the 
monochromatic luminosities (see also recently Sambruna et al. 2010). This trend results in 
a tight anti-correlation between the photon frequencies of the first peaks of the SEDs and 
the monochromatic luminosities. However, the frequencies of the first peaks are difficult to 



1 Two nearby GRBs, GRB 980425 and GRB 060218, are excluded from the sample to avoid the contami- 
nation caused by the radio emission from the associated supernovae. 
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be derived for the GRB afterglows at the current stage, because the SED modeling requires 
multi- wavelength follow-up observations from radio to hard X-rays (even 7-rays). To compare 
the GRB afterglows with the Blazar sequence, we use the radio-to-optical spectral slope a^o 
as a proxy of the frequency, both because q;ro is correlated with the frequency for the Blazars 
(e.g., Figure 8 in Fossati et al. 1998) and because o:ro can be easily derived for the GRB 
afterglows. The broad-band spectral slope a^o is estimated by following the definition 

«RO = ; - f l—Y- I 1 ) 

log(z/i/z/ 2 ) 

for both GRB afterglows and Blazars, where v x =5 GHz, A 2 = c/V 2 = 5500A (c is the light 
velocity). Similar as done in the above section, the spectral slopes a = —1/3 and a — 1.2 
is used to carry out /c-correction in the radio and optical band for the GRB afterglows, 
respectively (Sari et al. 1998). 

Figure 1 shows the estimated a RO plotted against the luminosities at 5 GHz. The Blazars 
from Fossati et al. (1998) are plotted by the red solid circles, and the GRB afterglows from 
Table 1 by the blue open ones. As expected from the Blazar sequence, there is a tight 
correlation between a^o and vL v for the Blazars. A spearman rank-order statistical test 
yields a correlation coefficient r s = 0.736 at a significance level P < 10 -4 where P is the 
probability of null correlation. The main conclusion drawn from the figure is that the radio- 
to-optical emission from GRB afterglows closely follow the Blazar sequence. The afterglows 
are clustered at the low luminosity end of the sequence only. The radio luminosities at 
5 GHz of the GRB afterglows cover a range from 10 40 to 10 42 erg s _1 , and a RO from 0.1 to 
0.5. According to the same statistical test, the correlation coefficient is slightly enhanced to 
r s = 0.798 with P < 10~ 4 when the two classes of the objects are merged into a single sample. 
As an additional illustration, Figure 2 plots the SEDs of the Blazar spectral sequence along 
with that of the afterglow of GRB 030329 observed at ~ 6day since the trigger. The average 
Blazar SEDs binned according to the radio luminosity were provided by Fossati et al. (1998) 
again. One can see directly from the figure that the SED of GRB 030329 follows the Blazar 
sequence at the low luminosity end. 

By combining the GRB afterglows and Blazars, an unweighted least-square fitting yields 
a relationship 

a RO = (0.27 ± 0.02) + (0.09 ± 0.01) log ( ^( 5GHz h (2 ) 

V10 4Z erg s V 

with a la dispersion of 0.09. 
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4. IMPLICATIONS 

The Blazar sequence was theoretically explained by Ghisellini et al. (1998). Only the 
very high energy electrons cool efficiently in a less powerful jet through synchrotron emission, 
which results in a break in the SED at high frequency. In contrast, an external seed photon 
injection (e.g., broad-line emission from AGNs) causes the powerful jet could efficiently cool 
for low energy electrons through the inverse Compton process. According to the unifying 
scheme, the low luminosity end of the vL v (5GH.z)-ano correlation is occupied by the high 
frequency-peaked BL Lac objects (HBLs). Without strong external radiation field, the jets 
in HBLs mainly cool through the synchrotron emission for low energy photons and weak 
SSC for very high energy photons, which predicts a very small Compton dominance, usually 
log(L c /L syn ) < 0. 

Our study indicates that the GRB radio/optical afterglows follow the Blazars in the 
^L [/ (5GHz)-q;ro diagram, which implies the afterglows share the similar radiation process as 
that occurring in the HBLs, that is mainly the synchrotron cooling without external photons. 
In fact, the weak inverse Compton components at X-ray band have been already identified 
in a few GRB afterglows with well multi- wavelength detections from radio to X-ray (e.g., 
Harrison et al. 2001; Frail et al. 2003b). The synchrotron model predicates that the peak 
photon frequency is f pea k oc -Br7p eak , where B is the magnetic field, Y the Lorentz boost 
factor, and 7 pca k the energy of the electrons emitting at the peaks of the SED. Assuming 
the recently updated theoretical Blazar sequence 7 pea k oc U' 1 (U oc L/(r 2 9 2 ) is the radiation 
energy density, where L is the luminosity, and r and 6j is the distance from center and half 
opening angel of the jet, respectively) for low luminosity Blazars (Ghisellini 2003) is still 
available for GRBs yields a relationship f pea k oc BTr A 6j / L 2 , which implies BTr 4 8j ~ Br A 8? 
(r ~ might be a constant for GRB late time afterglows taking into account of the tight 
Vpeak-vLv anti-correlation (see Figure 7 in Fossati et al. 1998). 

According to the widely accepted unification scheme for radio-loud sources, the weak 
Blazars are the beamed counterpart of FR I radio galaxies (Urry & Padovani 1995). The uni- 
fication scheme motivate us to suspect the radio emission of GRB afterglow is still detectable 
even though the strongly beamed GRB high energy emission does not point us. These after- 
glows without high energy prompt emission are called as orphan GRB afterglows. In fact, 
Bower et al. (2007) claimed that they identified a single candidate of radio orphan afterglow 
by examining the radio transient sources observed by VLA. 

The nature of optically "dark" bursts is still an open issue (e.g., Lamb et al. 2004). The 
dark bursts are now empirically defined as the events whose afterglow flux ratio of optical 
bands to soft X-ray at 3keV, /3 x, is less than 0.5 (Jakobsson et al. 2004). At the present, 
the debate on the optical darkness mainly focuses on two plausible scenarios: the afterglow 
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optical emission is intrinsic faint due to low efficiency (e.g., Jakobsson et al. 2005; Urata et 
al. 2007) or is heavily extincted by the dust (e.g., Zheng et al. 2009; Hashimoto et al. 20f 0). 

It is emphasized that none of the GRBs listed in Table I can be classified as a "dark" 
burst. Table 2 lists the radio and optical properties of three typical "dark" GRBs: GRB 970828, 
GRB 990506, and GRB 020819. For each burst, we list the GRB identification, measured 
redshift, observing frequency in radio, observation time elapsed since the trigger, observed 
flux density in radio, upper limit of R-band flux density, and the calculated monochromatic 
luminosity at the corresponding radio frequency. The data are compiled from literature as 
well. Along with the k- correct ion, 1) the flux densities (and also vL v ) at 5 GHz are trans- 
formed from that at 8.5 GHz by assuming a powerlaw spectral shape with slope a = — 1/3 
for GRB 990506 and GRB 020819; 2) the upper limits of R-band flux densities are trans- 
formed to V-band values by assuming a spectral slope a = 1.2 for all the three events. The 
three bursts are then over-plotted in Figure 1 by the green open triangles. The over-plotted 
vertical arrows show the calculated spectral slope a^o are underestimated because of the 
R-band upper limits. 

One can clearly see from Figure 1 that all the three "dark" bursts significantly deviate 
from the well established correlation (both optically bright GRB afterglows and Blazars) at 
a significant level larger than 2a due to their large ctcm- The deviation means these bursts 
are faint in optical bands with respect to not only X-ray emission, but also radio radiation. 

Briefly, there are two approaches that are widely adopted to interpret the optically 
"dark" burst puzzle: one is to reduce the value of f3ox by enhancing the X-ray emission, 
including long-lived reverse shock for X-ray (e.g., Dai 2004; Genet et al. 2007), late X-ray 
prompt emission (Ghisellini et al. 2007), and inverse Compton scattering (Pannaitescu 2008); 
an alternative approach is to suppress the optical radiation through heavy dust extinction 
or intrinsic faintness in optical bands. The larger «qr suggests that both X-ray and radio 
emission should be enhanced simultaneously for the optically "dark" bursts in the first 
approach. In contrast, the faint optical emission in the "dark" bursts could be plausibly 
explained by the heavy dust extinction (e.g., Djorgovski et al. 2001), because the radio 
emission is believed to be free of the dust extinction. The extinction scenario was recently 
supported by Zheng et al. (2009 and references therein) who found the natural hydrogen 
column density is systematically higher in optically "dark" bursts than optically bright 
events. The heavy extinction is also identified in the host galaxies of several "dark" GRBs 
(e.g, Levesque et al. 2010). 

In the context of the dust extinction scenario, we can roughly estimate the required 
extinction according to the a^o — vL v correlation. Taking into account of the definition 
of the slope q:ro oc \og(f Ul / f U2 ), we can express the GRB afterglow extinction Ay as a 
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following simple equation Ay = 2.5Aq;ro- Aoiro = «ro ~ ®ro * s the change of the slope due 
to the dust extinction, where a^o is the observed slope and a' R0 the expected one from the 
corresponding monochromatic luminosity according to Eq. (2) by assuming the o;ro — vL v 
relationship is universal for GRB afterglows. Note that this approach can only provide us the 
lower limit of the extinction. The la dispersion of the relationship gives a formal uncertainty 
of the estimated extinction 5Ay ~ 0.2 mag. The required extinctions are estimated for all 
the three "dark" GRBs and listed in Column (8) in Table 2. Generally speaking, the required 
extinctions A v > 0.5 — 0.6 mag, which is obviously larger than the averaged value (A v ) ~ 0.2 
obtained from a large sample of GRB optical afterglows (e.g., Kann et al. 2007). In fact, a 
typical dust extinction A v ~ 1 mag is identified for some optically "dark" GRBs (e.g., Perley 
et al. 2009). Recent near infrared observations indicate that the line-of-sight of afterglow of 
"dark" GRB 080325 shows heavy extinction with an amount Ay > 2 mag (Hashimoto et al. 
2010). 



5. SUMMARY 

We compile a sample of 18 GRB afterglows with multi-wavelength detections from 
radio to optical bands. Comparing the afterglows with Blazars shows that the afterglows are 
well consistent with the firmly established Blazar sequence (i.e., the vL v (5GHz)-o:ro The 
afterglows are only distributed at the low luminosity end of the sequence, which suggests 
that the GRB afterglows have similar emission mechanism with the high frequency-peaked 
BL Lac objects. The implication on "dark" GRBs is also discussed. 

The authors would like to thank the anonymous referee for his/her comments that im- 
prove the paper. This work was supported by the National Science Foundation of China 
(under grant 10803008) and by the National Basic Research Program of China (grant 
2009CB824800). 
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Fig. 1. — Radio-to-optical Broad band spectral slope «ro plotted against the monochromatic 
luminosity at 5 GHz at rest-frame. The Blazar sequence is plotted by the red solid circles. 
The 18 optically bright GRBs are over-plotted by the blue open circles, which are consis- 
tent with the sequence at the low luminosity end. The solid line presents the unweighted 
least-square fit to the data when the GRB afterglows and Blazars are merged. The la dis- 
persion is marked by the two dashed lines. In particular, the three green open triangles and 
corresponding vertical arrows mark the lower limits for the three typical optically "dark" 
GRBs (i.e., GRB 970828, GRB 990506, and GRB 020819). The "dark" GRBs deviate from 
the sequence by their large spectral slopes at a significance level > 2a. 
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Fig. 2. — The solid blue circles show the SEDs from radio to soft X-ray of GRB 030329 
afterglow detected at ~ 6 day since the trigger. The multi-wave length data are taken from 
5GHz, 8GHz, 15GHz, 23GHz, 43GHz: Berger et al. (2003), 90GHz: Kuno et al. 2005, 
250GHz: Sheth et al. 2003, 350GHz: Smith et al. 2005, V-band: Lipkin et al. 2004, 1 keV: 
Willingale et al. 2004. The other symbols plot the average SEDs of the Blazar sequence 
binned according to the radio luminosity, where the SEDs are taken from Fossati et al. 
(1998). 



Table 1. Radio/optical afterglow emission of the 18 optically bright GRBs 



GRB name 
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Table 2. Radio/optical afterglow emission of the three typical optically "dark" GRBs 



GRB name z At v F v Fb, vh v Ay Reference 

days GHz ^Jy ^Jy 10 40 ergs s^ 1 mag 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 



970828 
990506 
020819 



0.9578 5.6 4.86 99 < 0.5 
1.3 1.66 8.46 447 < 1.2 
0.41 1.5 8.46 315 < 4.4 



1.5 
1.3 
0.74 



> 0.63 Groot et al. 1998 

> 0.58 Taylor ct al. 2000 

> 0.50 Jakobsson ct al. 2005 



